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ABSTRACT 

Context. 1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high- 
energy (VHE, E > 100 GeV) y-rays due to its high X-ray and radio flux. 

Aims. Observations in the VHE y-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic 
particles present in the source and the radiation processes at work. Because of the distance of the source, the y-ray spectrum might provide further 
limits on the level of the Extragalactic Background Light (EBL). 

Methods. We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov 
telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are 
also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with 
the ATOM telescope are also used. 

Results. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is 
measured, corresponding to a significance of 7.8cr. The photon spectrum of the source is well described by a power law, with photon index of Tvhe 
= 3.45 ± 0.25 stat ± 0.20 syst . The integral flux above 200 GeV is (1.88 ± 0.20 stat ± 0.38 sysl ) xl0~ 12 cirr 2 s _1 . Observations with the Fermi-LAT in the 
first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 ± 0.2 stat ) x 10~ 9 ergcirT 2 s -1 , and a spectrum well described by 
a power-law function with a photon index The = L85 ± 0.18. SwiftfXKT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) x 
10~" ergcirT 2 s~', and a steep spectrum Tx = (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the 
synchrotron peak around 0. 1 ke V. 

Conclusions. Although the GeV-TeV observations do not provide better constraints on the EBL than previously obtained, they confirm a low 
density of the EBL, close to the lower limits from galaxy counts. The absorption-corrected HE and VHE y-ray spectra are both hard and have 
similar spectral indices (« 1.86), indicating no significant change of slope between the HE and VHE y-ray bands, and locating the y-ray peak 
in the SED above 1-2 TeV. As for other TeV BL Lac objects with the y-ray peak at such high energies and a large separation between the tw<p 
SED humps, this average broad-band SED represents a challenge for simple one-zone synchrotron self-Compton models, requiring a high Doppler 
factor and very low B -field. 
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1. Introduction 

The BL Lac object 1ES 0414 +009 was first detected with the HEAO 1 
satellite dGurskv et all 1 19781) in the energy range 0.2 keV-10 MeV, 
and i dentified in Einstein Observatory X-ray images dGiacconi et alj 
119791) . On the basis of the original X-ray, optical and ra dio observa- 
tions, it was identified as a probable BL Lac object by lUlmer et alj 
dl983l) . Located at a redshift of z = 0.287 dHalpern et al.1 Il99ll) . 
1ES 041 4+009 harbours a super-massive black hole of mass ~ 2 x 
10 9 M^ dFalomo et aljf2003l) . The hos t galaxy is classified as elliptical, 
with absolute mag nitude M R = -2 3.5 (Falomo et al. 2003). Polarization 
measurements by llmpev & Tapial dl988l) confirmed the classifications 
of this source as a BL Lac objec t. According to the classification scheme 
of iPadovani & Giommil dl995h . 1ES 0414+009 belongs to the class of 
high-frequency-peaked BL Lac (HBLs), objects with a synchrotron- 
emission peak located at UV/soft-X-ray frequencies, or equivalently, 
sources for which the X-ray emission is dominated by synchrotron ra- 
diation. 

There are several archival measurements of the spec- 
t ral properties o f 1ES 041 4+009 in X-rays tak en using ASCA 
Kubo et all 1 19981). RO SAT dBrinkmann et all 1 19951) . a nd BeppoSAX 



Beckmann et al. 



20021). Ob servations with EXOSAT dGiommi et al] 



990t ISambruna et alJI 1994h showed evidence of temporal variability 



and strong spectral variations, at times having an X-ray photon index 
of Tx < 2. During these states 1ES 0414+009 is characterized by 
a SED typical of "extreme" BL L ac objects (i.e. objects with the 
synchrotron peak above a few keV, ICostamante et al. [|2001|) . In the 
VHE y-ray domain, the data from the HEGRA experiment were used 
to derive an upper limit on the flux for thi s source correspondin g 
to 13.5 xl0~ 12 cm" 2 s _1 above 910 GeV dAharonian etafl 120001) . 
1ES 0414+009 was considered a good candidate for VHE emission by 
ICostamante & Ghisellini I d2002l) . on the basis of a high X-ray and radio 
flux, and its detection was yet more likely after blazar y-ray spect ra 
indicated a low intensity of the diffuse EBL (Aharonian et al. 2006a). 

Based on these VHE estimates, observations of 1ES 0414+009 with 
H.E.S.S. began in 2005 and continued until 2009, and significant ob- 
serving time was dedicated to its detection, because its high redshift 
made this source a potentially interesting candidate for EBL studies. 
In Sect. |2] the analysis of the 5-year data set on 1ES 0414+009 col- 
lected by the H.E.S.S. collaboration is reported. In the high-energy 
y-ray domain, emission from 1ES 0414+009 was first detected by 
the Large Area Tele scope (LAT) on the Fermi y-ray space telescope 
dAtwood et~alll2009l) . 1ES 0414+009 was detected with a significance 
of TS = 2Alog_£ > 25 in the firs t 1 1 months of the Fermi-LAS 1 oper- 
ation and is listed in the first year dAbdo et al. 2010a) and second year 
dAbdo et al.l 1201 ll) Fermi catalogues. In Sect. [3] the 20 month Fermi- 
LAT observation of 1ES 0414+009 is reported. The y-ray observations 
have been complemented with measurements in X-ray and UV, car- 
ried out using the SwiftP iKT and UVOT instruments dBurrows et alj 
1 20051 : iRoming et alj|2005l) . These, together with the o ptical measure- 
ments from the ATOM telescope dHauser et alj 120041) are presented. 
Constraints on the EBL are discussed in Sect. [4] and finally, the spectral 
energy distribution (SED) of 1ES 0414+009 is presented in Sect.[5] 

In this paper the luminosity distance of the source, located at a 
redshift z = 0.287 is computed with the standard ACDM cosmologi- 
cal model (H = 71 km/s/Mpc, C1 A = 0.73, Q. M = 0.27), yielding a 
luminosity distance D L = 1469.3 Mpc. 



2. VHE y-ray observations with H.E.S.S. 

H.E.S.S. is an array of four imaging atmospheric Cherenkov tele- 
scopes situated in the Khomas highland of Namibia (23° 16' 18" south, 
^"SO'OO" east), at an elevation of 1800 metres above sea level (see 
lAharonian et al. I2006bl) . The H.E.S.S. array is sensitive to y-rays above 
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Fig. 1. The distribution of the squared angular distance (6 2 ) 
for the ON-source events (crosses) and the OFF-source events 
(shaded). The dashed blue line represents a PSF profile for a 
point-like source estimated by Monte Carlo simulations. 



~ 100 GeV, and commonly achieves an event-by-event angular resolu- 
tion of <0. 1° and a relative energy resolution of ~15%. 

The BL Lac 1ES 0414+009 was observed every year between 
2005 and 2009. After run-quality selection the data set comprises 73.7 
hours (live time) of observations, of which 67.5 hours were taken 
with four telescopes. Observations have been carried out at zenith an- 
gles of 22° to 41° with a mean value of 26°. The pointing offset 
from 1ES 0414+009 was 0.5°. The dat a were analysed with the Model 
Analysis (de Naurois & Rolland 2009), in which shower images of all 
triggered telescopes are compared to a precalculated model by means 
of a log-likelihood minimis ation. The analys is has been cross-checked 
with a multivariate analysis dOhm et al.l2009t) using an independent cal- 
ibration, which yielded consistent results. 

Two sets of cuts were used. For the source-detection analysis stan- 
dard cuts were used, including a minimum image charge of 60 photo- 
electrons (energy threshold E th = 220 GeV, defined here as the energy at 
which the effective area falls to 10% of its maximum value). Loose cuts, 
with a smaller charge cut of 40 photoelectrons, result in a decreased 
signal-to-background ratio but have the advantage of higher photon 
statistics and a lower threshold of 150 GeV. Since 1ES 0414+009 is a 
very faint and distant source, loose cuts were used for the spectral anal- 
ysis and light curve generation Q Therefore, in this text all the fluxes 
and light curves shown are derived with loose cuts. 

A significant excess of 224 events (7.8<r) from the direction of 
1ES 0414+009 was found in the total data set using standard cuts. A 
year-by-year summary of the observations and of the corresponding re- 
sults is shown in Table [T] The background subtracted distribution of 
the square of the angular difference between the reconstructed shower 
position and the nominal source position (6 2 ) is shown in Fig. [T] The 
Reflected-Region method dBerge et al.l2007l ) was used for the definition 
of the ON-source and OFF-source data regions. The 2 distribution is 
consistent with that of a point-like source, as shown in Fig.[T] 

The fit of the uncorrected excess map with a two-dimensional 
point-like function convolved with the instrument point spread function 
PSF (68% containment radius of 0.064° for standard cuts) also indi- 
cates that the source is consistent with being point-like and is located at 
(aj2ooo=4 h 16 m 52.96 s ± 0.10^, ± 0.10 s , < 5 J2 ooo=l 5'20.4" ± 15^). This 
position is consistent with t he nominal positio n (ffj 2 ooo=4 h 16 m 52.8 s , 
<5j20oo = 1°5'24") reported bv lUlmer et aLldl983l) . 

Fig. [2] shows the differential energy spectrum of the VHE y-ray 
emission above the energy threshold of ~ 150 GeV using the ~530 ex- 
cess events obtained with loose cuts. The spectrum is obtai ned using a 
forwar d-folding maximum-likelihood method, described bv lPiron et al.1 
d200ll) . The spectrum is compatible with a power-law distribution 



1 In this case the loose cuts almost double the number of excess 
events resulting in a poorer background rejection. 
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Table 1. Summary of H.E.S.S. observations of 1ES 0414+009 for each year. The MJD of the first and last night of the observation 
in each epoch between 2005 and 2009, the live-time, the number of ON- and OFF- source events, the excess and its significance 
are reported. These events are obtained with standard cuts and for all the observation the normalization factor between the OFF and 
ON exposures is 12.6. The last three columns list the integral flux above 200 GeV, the x 2 , number of degrees of freedom (ndf) and 
the x 1 probability P(x 2 ) for a constant fit to the light curves binned by night within each year (first five rows), or yearly (last row) 
within the total observing period. 



Epoch 


MJD starl 


MJD stop 


Live-time[h] 


ON 


OFF 


Excess 


SigM 


F(E >200 GeV) ( *» 
10- 12 cm" 2 s" 1 


y/ndf ( "» 


Pix 1 ) l " } 


2005 


53646.06 


53706.97 


15.2 


162 


1474 


44.9 


3.8 


2.16+0.41 


20.5/16 


0.20 


2006 


54024.03 


54120.82 


11.4 


189 


1834 


43.3 


3.3 


1.19+0.51 


10.5/12 


0.57 


2007 


54406.98 


54475.84 


12.9 


150 


1541 


28.1 


2.4 


1.73+0.46 


13.9/11 


0.24 


2008 


54736.09 


54750.10 


9.1 


117 


1160 


25.1 


2.4 


1.37+0.62 


2.5/8 


0.96 


2009 


55067.14 


55158.97 


25.2 


289 


2602 


82.8 


5.2 


1.86+0.41 


31.0/28 


0.32 


Total 






73.7 


907 


8611 


224.3 


7.8 


1.84+0.20 


2.66/4 


0.62 



Derived with loose cuts and with the spectral index fixed to the global value Tvhe = 3.45. 
Night-by-night and year-by-year light curves are derived with loose cuts. 



(&NI&E ~ E- r ) with a photon index of r VH E = 3.45 ± 0.25 stat ± 0.20 syst ; 
the differential flux at the decorrelation energy of 305 GeV is (5.70 ± 
0.62 stat ) xl0~ 12 cnr 2 s-' TeV~'. In Fig[2]the confidence band of the 
power-law fit for 1ES 0414+009 is illustrated, clearly steep in appear- 
ance and affected by EBL absorption. The spectral points (with lcr sta- 
tistical errors) in the upper panel in Fig.|2]are derived from the residuals 
in the different energy bins shown in the bottom panel in the same fig- 
ure. The integral flux above 200 GeV is F = (1.88 ± 0.20 stat ± 0.38 syst ) 
xl0~ 12 cirr 2 s -1 and it correspon ds to ~ 0.6% of the Crab Nebula flux 
above the same energy threshold ( Aharonian et al. 20d6bT). 

The light curves for E > 200 GeV with different time samplings 
have been derived with loose cuts, assuming the spectral index r V HE = 
3.45 obtained by the maximum-likelihood fit (see the top panel in 
Fig. f3). No evidence of variability is found on different timescales. 
Accounting for both statistical and systematic (20%) uncertainties on 
the flux points, the fit of the night-by-night light curve with a constant 
value yields x 2 = 81.1 (for 79 degrees of freedom) with a correspond- 
ing probability P=0.41 and a 99% confidence level upp er limit on the 
fractio nal variability F var < 0.40 with F var defined as in IVau" ghan et alj 
d2003l) (calculated with the method described in iFeldman & Cousins! 
1998). The fit of the monthly light curve yields a x 2 = 18.73 for 12 
degrees of freedom, with a corresponding probability of 0.1 and upper 
limit F var < 0.76. The same conclusion is found on a yearly timescale, 
with a probability of a constant flux P=0.62 and 99% confidence level 
upper limit F var < 0.32. 

Table [T] shows the measured integral flux for each year the source 
was observed. Also reported are the x 2 values for a fit of a constant 
when the flux is binned by nights within each year and by years within 
the total observation. As the x 2 probability for each fit is greater equal 
than 0.2 there is no evidence for variability in this data set. 

3. Multi-wavelength observations 

3. 1. Gamma-ray observations with the Fermi LAT 

The LAT is a pair-conversion telescope, and is the primary instrument of 
the Fermi satellite la unched in June 2008 . The Fermi-hAT instrument, 
described in detail in lAtwood et a l. ( 2009), detects y rays with energies 
between 20 MeV and >300 GeV. The bulk of the Fermi observational 
program is dedicated to a survey, in which the full y-ray sky is covered 
every 3 hours. 

The spectrum of 1ES 0414+009 in the 100 MeV to 100 GeV energy 
ban43 was derived from LAT data taken between August 4, 2008 and 
May 17, 2010 (651 days). To minimize background contamination, only 
events with a high probability of being y rays, known as "DIFFUSE" 

2 The analysis has been restricted to the energy range over which the 
response of the instrument is best understood. 
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Fig. 2. Differential energy spectrum of 1ES 0414+009. The 
shaded band corresponds to the range of the power-law fit 
(68% confidence level), taking into account statistical errors. 
The lower panel shows the residuals of the fit; the difference 
in each energy bin between the expected and observed number 
of excess events, normalized to the uncertainty on the latter. The 
6 8% confidence level upper l imit is calculated with the method 
of lFeldman & Cousins! (fl998b . 



events, having an angle of less than 105° to the Zenith, were retained. 
The spect rum was evaluated using an unbinned maximum likelihood 
approach (Mattox et al. 1996) in which the events originating within a 
region of interst (ROI) of 10° around 1ES 0414+009 were modelled us- 
ing diffuse Galactic and isotropic components and 17 point sources with 
power-law spectra. The diffuse components are tabulated in the standard 
files glLiem_v02.fit and isotropicjem_v02.txt, respectively; the overall 
normalization of each component was allowed to vary in the fit. The 
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point sources include a test source at the l ocation of 1ES 041 4+009, 
all sources from the first Fermi catalogue dAbdo et ail 1201 0j) within 
12° of 1ES 0414+009 and three additional sources found in a dedi- 
cated search for additional y-ray emission in the ROI. These additional 
sources are weak, with detection significances between TS = 24 and 
TS = 42 (approximately Act to 5.6cr), and all lie relatively far from 
1ES 0414+009 , the closest being separated by > 5.213 For the brighter 
background sources the power-law spectral assumption was tested by 
evaluating a log parabolic model, which was found not to significantly 
improve the fit. In total the model had 30 free parameter^ The maxi- 
mum likelihood optimization was performed using the standard Fermi 
ScienceTools package (v9iT9pO), using the P6_V3_DIFFUSE instru- 
ment response functions (IRFs). 

The null hypothesis, that no source is present at the location of 
1ES 0414+009, is rejected with TS = 67.7 for 2 degrees of freedom, 
or approximately 7.9cr. The best-fit power law has an integral flux of 
(4.3 ± 2.2 stat ^ ) x 10~ 9 cm" 2 s" 1 (100 MeV to 100 GeV) and an 
index of The = 1-85 ± 0.18. The differential flux at the decorrelation 
energy of 3.6 GeV is (5.0 ± 1.1) x 10~ 14 cnr 2 s~> MeV-'. The sys- 
tematic errors on th e index, estimated using the bracketing IRE method 
dAbdo et al. 2010b), in which the IRFs are modified to give the greatest 
change in the index of the best-fit power law, are THE.syst =1q 

The most energetic photon detected within a region of 0.5° around 
the source was reconstructed with an energy of E = 19.7 GeV. The 
energy resolution is less than 10% for photons with E = 20 GeV. 

Variability in the flux of y rays on monthly timescales was eval- 
uated by producing a l ight cu rve with 32.5 day bins (second panel of 
Fig. [3] see lAbdo et alJd2010al) for details of methodology). Comparing 
the likelihood of a model in which the flux in each time bin is free, 
to one where it is assumed constant gives a difference of TSvar = 
2AlogX = 19.2 for 19 extra degrees of freedom (DOF). The prob- 
ability of getting a valu e of TS V ar > 19.2 by chance is P = 0.44, 
assuming the theorem o f lWilksl dl938h is applicable, and we conclude 
that there is no evidence for variability on monthly timescales. Since 
there is a relatively large number of upper limits in the 32.5-day light 
curve, we rebinned the data to produce a 130-day light curve which is 
also shown in Fig. [3] (blue points). Applying the same variability test 
gives TS VAR = 1.44 for 4 DOF, with a chance probability of P = 0.84, 
again consistent with no variability. 

3.2. X-ray observations with the Swift/XRT 

Archival Swift observations within the timespan of the overall H.E.S.S. 
observing campaign were analysed. Swift performed three snapshot ob- 
servations in 2006, on October 21 (1138 s), 22 (226 s) and 23 (553 s), 
and three in 2008, two on January 29 (1008 and 993 s) and one on 
February 4 (2227 s). 

The XRT data were processed with standard procedures using the 
FTOOLS task XRTPIPELINE (version 0.12.4). Only data taken in pho- 
ton counting (PC) mode were considered, given the low rate of the 
source (< 0.5 counts/s in the 0.3-10 keV range) and because the ex- 
posures in Window Timing mode were extremely short. Source events 
were extracted in the 0.3-10 keV range within a circle with a radius 
of 20 pixels (47"), while background events were extracted from both 
circular and annular regions around the source, free of other sources. 

The spectra were fitted with a single power-law model with Galactic 
absorption N u fixed at 9. 15 x 10 20 cm" 2 , as obtained by dedicated 21 cm 
observations dElvis et akll 19991) . No significant flux or spectral variabil- 
ity was found within the 2006 or Jan 2008 data sets (the fit with a con- 
stant gives ^ 2 /ndf of 4.9/10 and 10.5/7, respectively, with a probability 

3 These additional sources are not associated with any point sources 
in 2FGL, however one is approximately two degrees from a 2FGL 
source marked as confused. It seems likely that they are the result of 
deficits in the galactic diffuse model that were eliminated in the im- 
proved model use in 2FGL. 

4 The spectral parameters for point s ources formally out side the 10° 
ROI were fixed at their catalogue values dAbdo et alj2010al) . They were 
included in the model to account for the effects of the relatively large 
PSF at low energies. 



of 0.90 and 0.16). The corresponding exposures were thus summed, ob- 
taining a total of three data sets for further analysis: the first covering 
the period October 21-23, 2006 (1916 s), the second for January 29 
(2001 s), and the third for February 4, 2008 (2227 s). 

To test the possible presence of spectral curvature, broken power- 
law and log-parabolic models were compared to the simple power-law 
hypothesis. The three spectra are well fitted by a single power law 
within the available statistics, with steep indices around 2.2-2.3 (see 
Table [2). Log-parabolic models do not significantly improve the x 1 
(F-test~85-95%). Variations among the different epochs are modest: 
both flux and spectra are similar between October 2006 and January 
2008, while on Feb 4 the flux increases by ~15%, but consistent within 
1 sigma with other two epochs, with a slightly harder spectrum by 
Ar 0.15 ± 0.10. The light curve for the XRT observations is shown in 
Fig.[3] 

The Swift Burst Alert Telescope (BAT. lBarthelmv et al. 1120051) data 
from 200 4 to 2010 were also analysed, following the procedure de- 
scribed in lAielloetaT] d2008h . The' source presented only a marginal 
signal of 2.5cr in the 14-30 keV band. The corresponding 99% upper 
limit on the flux is reported in the SED in Fig. [4] for reference. 

3.3. Optical and UV observations with the Swift/UVOT 

The UVOT instrument dRoming et alll2005h took exposures in differ- 
ent filters during the XRT pointings, namely V (547 nm), B (439 nm), 
and U (347 nm) in October 2006, UVM2 (225 nm) in January 2008 
and U in February 2008 (the observations are summarised in Table|2j. 
Photometry of the source was performed using the UVOT software in 
the HEAsoft 6.9 package. Counts were extracted from the standard 
aperture of 5" radius for all single exposures and all filters, while the 
background was carefully estimated from different positions more than 
27" away from the source. Count rates we re then converted to fluxes 
using the standard photometric zero points dPoole et al.ll2008h . 

The fluxes were d e-reddened fo r Galac tic absorption using the 
extinction curve from ICardelli et alj d!989l) with the updates from 
IO'Donnellldl994l) . assuming R w [= A V /£(B- V)] = 3.1. This is the av- 
erage value for the Galactic diffuse interstellar medium. For the line of 
sight of 1ES 0414+009, a value of A B = 0.507 was used dSchlegel et aTl 
119981) , wifh£(B - V) = 0.117 mag. 

The host galaxy of 1ES 0414+009 is elliptical, of total R magni- 
tude 17.49 and half-light rad ius 4.7" ± 0.5" (from HST snapshot ob- 
servations, Scarpa et al. 2000). To isolate the flux of the active core the 
UVOT fluxes were corrected for the contribution of the host galaxy. The 
wavelength-depen dent correction w as determined using a template for 
elliptical galaxies dSilva et al.ll 19981) . rescaled to the host-galaxy flux in 
the R band, and accounts for the given apertures. 

As in the X-ray band, no significant variability was found in the UV 
band in October 2006, nor in January 2008 (see fourth panel in Fig. [3] 
where fluxes are reported without extinction correction). The individ- 
ual frames were thus co-added, and the results are given in Table [2] A 
decrease by 0.5 mag can be noted between 2006 and February 2008, 
in the U filter (the only filter in common). This behaviour contradicts 
what observed in the X-ray band, where the corresponding XRT spec- 
trum showed an increased flux and harder spectrum. This is indicative 
of a shift of the SED peak towards higher energies, which is confirmed 
by the simultaneous UVOT-XRT spectral fitting (see Sect. 14. lb . 

3.4. Optical observations with ATOM 

The 75 cm optical telescope of the H.E.S.S. collaboration, ATOM, is lo- 
cated next to the Cherenkov telescopes at the H.E.S.S. site in Namibia. 
ATOM has been monitoring 1ES 0414+009 since Nov 2006 with a sam- 
pling of »2 observations per week during its visibility from July to 
March. 

During the H.E.S.S. observations of 1ES 0414+009, the sampling 
was increased to obtain truly simultaneous data. The observations 
of this source were mostly carried out in B and R band (440 nm 
and 640 nm). On some of the nights, observations in V and I bands 
(550 nm and 790 nm) were also performed. The data analysis con- 
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Table 2. Spectral parameters and fluxes of the Swift data in the three epochs. 



Quantity 


21-23 Oct 2006 


29 Jan 2008 


4 Feb 2008 




Units 




54029.04-54031.04 


54494.23 


54500.21 




MJD 


Single power-law fit, XRT data 


r 


2.30 ± 0.08 


2.33 + 0.08 


2.18 + 0.07 






^2-lOkeV 


8.41 ± 0.70 


8.92 ± 0.70 


10.3 ± 0.7 


io- 


■12 -2 -1 

erg cm s 


*>df 


1.06/42 


1.00/22 


0.80/43 






UVOT filters, observed magnitudes 1 "' 


V 


16.56 + 0.03 










B 


16.97 ± 0.03 










U 


16.06 ± 0.03 




16.51+0.03 






UVM2 




16.56 + 0.03 










Log-parabola 1 b) fit, 


combined UVOT + XRT data 






UVOT filters 


B,U 


UVM2 


U 






TlkeV 


2.28 ± 0.07 


2.29 ± 0.08 


2.16 + 0.06 






b 


0.12 + 0.03 


0.17 + 0.03 


0.10 + 0.03 




curvature 


peak 


0.071 +0.007 


0.13+0.02 


0.17 + 0.06 




keV 


^2-10keV 


7.83 + 0.70 


8.39 + 0.70 


9.96 ± 0.70 


10' 


■12 -2 -1 

erg cm - s 


*>df 


1.01/43 


0.84/23 


0.79/44 







Errors are reported at 90% confidence level for 1 parameter (&x~ = 2.71). 
Statistical errors from uvotsource analysis. 
<b > Log-parabola model is defined as: W/dE k £- r+ft l0 s< £ > 

(c) Value of t he peak energy in the SEP, obtained by fitting the same data with a different parameterization of the log-parabolic model, for which the free parameters are E^ik and b 
(details in Tramacere et al. 2007). 



sists of debiassing, flat fielding and photometry using SExtractor 
dBertin & Arnou ts 1996) a nd was perfo rmed using the ATOM-specific 
pipeline processing dHauser et al.ll2004h . The fluxes were determined 
using a 4 radius aperture and the flux scale wa s fixed on 7 nearby stars 
from the USNO catalogue dFiorucci et al.lll998h . 

The night-by-night light curve of the ATOM observations with the 
R-filter is shown in the bottom panel of Fig. [3] where no Galactic ex- 
tinction correction was considered. There is no strong variability in the 
light curve (the amplitude variations are less than 50% overall) but a 
long-term trend can be seen, with a minimum in 2008 with respect to 
the other years. No change in (B-R) colour was observed. 

4. Combined spectral fits 

4. 1 . UVOT and XRT spectral fits 

To study the connection between the UV and X-ray components of the 
synchrotron emission, and to better estimate the position of its peak, a 
broad-band fit of the simultaneous UVOT and XRT data was performed 
in Xspec, using the UVOT response matrices in the calibration database. 
Only the B and UV filters were considered, to minimize the contami- 
nation from the host galaxy. In Xspec, Galactic reddening and X-ray 
abso rption are taken i nto account using the models redden (updated 
with lO'Donnell 119941 ) and wabs for the two bands, respectively. 

The UV-to-X-ray broadband spectra in all three epochs are well 
represented by a single log-parabolic component (Table |2j, with val- 
ues of the curvature which are typical for HBLs: around 0.1-0.2 
dTrama cere et al. 2007[). The resulting estimate for the energy of the 
synchrotron peak is around 0.1 keV, overall, with a slight shift (at ~ 3cr 
level) of the peak position from 0.07 to 0. 17 keV between the 2006 and 
2008 epochs. The corresponding flux at the synchrotron peak is steady 
around (1.25 ± 0.10) x 10"" ergcirr 2 s~' between 2006 and Jan 2008, 
decreasing slightly to ( 1 .0 1 ± 0. 1 1 ) x 1 0~ 1 1 erg cm" 2 s" 1 on Feb 4, 2008 . 

4.2. Fermi and H.E.S.S. spectral fits 

By combining the Fermi and H.E.S.S. spectra it is possible to study the 
broadband properties of the y-ray SED peak, once the effects of y - y 
interactions with photons of the diffuse EBL are taken into account. At z 
= 0.287, the intrinsic spectrum calculated from the VHE observations is 
strongly sensitive to the EBL density. Using the high EBL density given 



by dire ct estimates at 1-3 //m (~20 nW/m 2 sr, see e.g. Aharonian et al. I 
I2006al) , the reconstructed spectrum is unrealistically hard (r illt < 0.3), 
but it softens to Ti nt > 1.5 for lower EBL densities. 

Considered independently of other spectral information, therefore, 
the H.E.S.S. spectrum of 1ES 0414+009 further corroborates the up- 
per limits on the EBL level previously derive d from other HBLs 
d Aharo nian et al. 2006a; Ahar onian et al.1 12007 al ibi), indicating a low 
density around ~l-3 yum close to the lower bound given by galaxy 
counts dAharonian et al. 2006a). However, the measured spectrum of 
1ES 0414+009 does not improve those upper limits, nor their statistical 
uncertainty. 

For example, comparing with 1E S 1101-232 (z = 0.186) a nd per- 
forming the same analysis described in Aharonia n et al. I d2006al) . where 
a scaled EBL spectral shape was used, a higher EBL scaling factor P is 
required in this case as compared to 1ES 1101-232 (P0.58 vs P0.55) 
in order to give the same estimated intrinsic spectral index. The higher 
redshift (z = 0.287 vs z = 0.186) provides a longer path-length over 
which the effect of EBL absorption is integrated, thus giving a larger 
measurable effect of absorption of the intrinsic spectrum (i.e. a smaller 
difference in EBL flux AF EBL is sufficient to yield the same observed 
Ar). This longer leverage, however, is not enough to compensate for 
the larger statistical error on the spectral slope (~0.35 vs ~0. 15). 

Therefore, the approa ch adopted here is to fix the EBL to the 
best current estimates (e.g. Franceschini et al. 2008; Dominsuez et al. 
1201 lh . and to investigate the resulting S ED properties. The EBL model 
assumed here ((Franceschini et al. 2008]) is close to the lower limits from 
galaxy counts and comp atible with the limits from VHE observations 
dAharonian et al. l2006al) . With this model, both the Fermi and H.E.S.S. 
absorption-corrected spectra are harder than 2, indicating a high-energy 
peak above few TeV Both spectra have the same index, within sta- 
tistical errors (The - 1-85 ± 0.18 vs r V HE,int — L87 ± 0.35), indicat- 
ing that the spectral properties do not change significantly between the 
HE and VHE bands. The typical uncertainty on the VHE spectral in- 
dex due to the residual EBL uncertainty can be estimated to be in the 
range [+0.15,0.3], scaling the Franceschini template between the level 
of galaxy counts and the current upper limits. A single power-law model 
provides a good representation of the combined Fermi and H.E.S.S. 
spectra from 1 GeV to 2 TeV, with index 1.86 ± 0.06 and x 1 jndf = 
4.7/7. 

Because the H.E.S.S. observations overlap only partially in time 
with the Fermi data, variations in the flux could bias the combined spec- 
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Fig. 3. Overall light curve of the multi-wavelength observations of 1ES 0414+009 covering the period from 2005 to 2010. Night- 
wise (shaded) and year-wise (dark) light curves of 1ES 0414+009 by H.E.S.S. (a). The SwiftfKRT (c), UVOT (d), and ATOM (e) 
observations have a night-wise binning, while the Fermi light curve (b) has 32.5 day binning with 95% upper limits (black open 
circles) and 130 day binning (blue filled circles). Only statistical uncertainties of the flux points are shown. 



tral fitting of the average spectra towards harder or softer values. To 
quantify this, the Fermi average spectrum during the H.E.S.S. observ- 
ing epochs (~3.6 months) was compared to the total (20 months). 

The differences are negligible with respect to the statistical errors 
(<7%), and the confidence contours in the flux-index plane of the total 
data set are fully enclosed in the respective contours of the H.E.S.S. 
epoch. Thus, it is possible to conclude that no correction between the 
Fermi and H.E.S.S. spectra is necessary. 



5. Discussion 

The overall SED is shown in Fig. [4] Unfortunately, the X-ray point- 
ings are not strictly simultaneous with the H.E.S.S. or Fermi observing 
epochs, (see Fig. f3}. However, given the lack of strong variability in 
all the observed bands, it can be reasonably assumed that the ATOM, 
Swift, H.E.S.S. and Fermi data provide together a good representation 
of the average SED properties of the source over the time span of the 
observations, 

These observations confirm 1ES 0414+009 as a typical HBL, with 
a synchrotron peak around 0.1 keV and where the X-ray band is domi- 
nated by the synchrotron emission of high-energy electrons. The higher- 
energy component, which is commonly interpreted as inverse Compton 
(IC) emission from the same population of electrons responsible for 
the synchrotron emission, was largely unexplored before the obser- 
vations reported here. The combined data indicate that a single, hard 
power law fits the intrinsic spectrum over the combined HE-VHE y- 



ray ranges, constraining the peak of the IC emission to lie at energies 
greater than -1-2 TeV. Thus, 1ES 0414+009 seems to belong to the 
subclass of HBLs characterized b y Compton peak energies in the multi- 
TeV range, like 1ES .1101- 232 dAharonian et al J l2007bl). IBS 0347- 
121 jAharonianetal|[2007ah . 1ES 0229+ 200 dAharonian et alj|2007ct 
iTavecchio et all2009h . or 1ES 1218+304 dAcciari et al.ll2009l) . 

Such high energies for the IC peak are difficult to explain with a 
standard one-zone SSC model, due to the decrease in the scattering 
efficiency in the Klein-Nishina regime and the decrease in the energy 
density of the seed photons available for scatterings in the Thomson 
regime. Both these effects tend to steepen the emitted y-ray spectrum 
at VHE. In principle it would be possible to avoid this steepening by 
assuming that the cooling of the electrons occurs only by IC scattering 
in the deep Klein-Nishina regi me (for example, off a black-body-like 
spectrum of target photons, see lModerski et al]|2005l) . However, the re- 
sulting hardening of the electron distribution at high energies (because 
high-energy electrons would cool slower than those at lower energies) 
would necessarily imply a stron g hardening of the sync hrotron spectrum 
towards hard X-ray energies dModerski et al ] 120051) . The Swift/XKT 
and BAT data do not seem to corroborate this hypothesis, though such 
an enhancement could in principle be located at even higher energies 
(> 100 keV). 

Therefore, to explain an SED such as that of 1ES 0414+009 (Fig.|4} 
in the context of a one-zone SSC model, one is drawn to consider a set 
of parameters for which t he IC peak is obtained i n the Thomson regime 
(jpeakhv' ^ m e c 2 /4) (see lAharonian et aD Ilxil iBoettcher et aT]|2008l. 
for alternative scenarios). The branch of the electron distribution corre- 
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sponding to the GeV-TeV band would then correspond to the optical- 
UV band of the synchrotron spectrum (i.e. electron energies y < y pi ,ak)- 

Adopting hv s =0.1 keV and hv c = 2 TeV as the energies of the 
synchrotron and IC peaks respectively, the Lorentz factor of the elec- 
trons at the peak would be y p = (3v c /4v s ) 1/2 =s 1.2 x 10 5 . The lu- 
minosity at these two energies can be estimated as L s ~ 10 45 5 and 
L c ~ 10 44 - 9 erg/s, respectively, from the log-parabolic fit to the syn- 
chrotron peak and the power-law fit to the y-ray spectrum. Performing 
an analysis of the independent SSC constraints in the (log B - log 6) 
plane (see e.g. iTavecchio et afll 19981 : iGuv et all 2000). with these pa- 
rameters, a solution in the Thomson regime can be found for 5 > 50 and 
B < 0.0 1 G, but only with a rather large emitting region (R ~ 10 17 cm, vs 
Ri- ~ 6x 10 14 cm, the Schwartzschild radius of the putative Black Hole). 
If this is the case, no variability faster than 1 day is expected, though 
this is below the sensitivity of current experiments at this source's flux 
level. Assuming instead a more typical size of the emitting region for 
HBLs, around R ~ 1 - 2 x 10 16 cm, rather extreme parameters would 
be required: 6 > 200 and B < 0.002 G. 

However, in both cases the synchrotron cooling is very slow, and 
the cooling time of the electron s at the synchrotron pea k is much longer 
than the escape time (see e.g. ITavecchio et alj[l998t considering the 
electron escape velocity in the range l-l/3c). The break in the elec- 
tron spectrum corresponding to the SED peak, therefore, cannot be 
explained simply by radiative cooling (a "cooling break") or by an 
equilibrium between cooling and escape from the source, as typically 
assumed for BL Lac objects with much lower Com pton-peak ener- 
gies like PKS 2155-304, Mkn 421 or Mkn 50 1 (e.g. ITavecchio et alj 
Il998l : iGuv et aU200Cl : iKrawczvnski et aU2002l : iGhisellini et aU2002h . 
Furthermore, the slow cooling requires the kinetic energy of the jet to 
be increased significantly in order to compensate for the low efficiency 
of the radiative emission. It should be noted that the Fermi and H.E.S.S. 
spectra provide only a lower limit on the energy and luminosity of the 
IC peak; if the true source values are actually higher, even more extreme 
parameters would be required. The problems just described are a gen- 
eral issue for the class of TeV BL Lac objects characterized by a hard 
TeV spectrum, an IC peak at very high energies and a large separation 
between synchrotron and IC peak frequencies. It might indicate either a 
fundamental difference in the physical jet conditions with respect to the 
other HBLs, or the need for an emission mechanism different from that 
of homogeneous SSC models. 

As a final caveat, however, it should be noted that the statisti- 
cal uncertainty on the y-ray spectra is rather large, especially at VHE 
(TvHE,int = 1.87 ± 0.35). If the true VHE spectrum were > lcr steeper 
than the H.E.S.S. measurement (i.e. r V HE,int = 2.2), the IC component 
would have a peak energy around 100-200 GeV, and a slightly lower 
luminosity. This would make the 1ES 0414+009 SED more similar to a 
typical HBL and would be easier to explain with standard SSC param- 
eters. 



6. Conclusion 

The H.E.S.S. array of Cherenkov telescopes has been used to detect sig- 
nificant VHE y-ray emission in the direction of 1ES 0414+009, a distant 
BL Lac objects of known redshift. The overall flux above 200 GeV from 
2005 to 2009 shows no variability on any time-scale and is on average 
very low (around 0.6% Crab Nebula flux), making 1ES 0414+009 one 
of the faintest extra-galactic sources detected in the TeV domain. This 
source was detected with the Fermi LAT in the first 20 months of its 
operation between 2008 and 2010 and was very faint in the HE domain 
too. 

The VHE spectrum measured using H.E.S.S. is consistent with the 
current limit of the EBL and confirms the low level of EBL in the 
few fjm range as derived from galaxy counts. The HE and VHE spec- 
tra (absorption corrected with an EBL model close to the lower limits) 
show a best-fit power law with an index harder than 2, indicating that 
1ES 0414+009 can be classified as a hard-TeV BL Lac object. The over- 
all SED of this source is averaged over five years with no strict simul- 
taneity between the Swift and H.E.S.S. observations. With this caveat, 
the properties of the SED -in particular an IC peak energy above 1- 



2 TeV- are difficult to explain in the framework of a pure one-zone SSC 
model, unless using unusual values for the main parameters. 

Future observations of 1ES 0414+009 are strongly motivated, in 
particular with the next generation of Cherenkov telescope arrays such 
as the Cherenkov Telescope Array (CTA). With a tenfold increase in 
sensitivity with respect to the current Cherenkov experiments, CTA will 
allow better measurements of the VHE spectrum of faint sources, espe- 
cially at high energies, provide better constraints of the IC peak, and 
will help to address the physics issues posed by the hard-TeV BL Lac 
objects. 
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Fig. 4. Average SED of 1ES 0414+009 from the observations reported in Fig. [3]in the epoch 2005-2009. H.E.S.S.: black filled/open 
circles with/without EBL correction, respectively. Fermi: red full circles correspond to detections with TS > 9, red downward- 
pointing arrows to 95% upper limits. Swift/BAT in 5 yrs: magenta 99% upper limit. Swift/XRT&UVOT: green open squares (Oct 
2006), blue open triangles (Jan 2008), and cyan crosses (Feb 2008). ATOM: black triangle and blue inverted triangles corresponding 
to two different states, lower flux level in 2008 and higher in 2006, respectively, in the R and I bands. Gray points and butterflies are 
a collection of archival data from lCostamante & Ghisellini I (120021) and references therein. 
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